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* „ " spacing between the subcarriers without a large degradation in 

system performance [I J. 

There have been several papers on (he subject of synchro- 
nization for OFDM in recent years. Moose gives the maximum 
likelihood estimator for the carrier frequency offset which is 
calculated in the frequency domain after taking the FFT [2J. 
He assumes that the symbol timing is known, so he just. has to 
find the carrier frequency offset. The limn of the acquisition 
range for the carrier frequency offset is ±1/2 the subcarrier 
spacing. He also describes how to increase this range by using 
shorter training symbols to find the carrier frequency offset. 
For example shortening the training symbols by a factor of 
two would double the range of carrier frequency acquisition. 
This approach will work to a poirit. but the estimates get 
worse as the symbols get shorter because there are fewer 
samples over which to average, and the training symbols need 
to be kept longer than the guard interval so that the channel 
impulse response does not cause distortion when estimating 
the frequency offset. 

Nogami and Nagashima [4] present algorithms to find the 
carrier frequency offset and sampling rate offset. Hjey use 
a null symbol where nothing is transmitted for one symbol 
period so that the drop in received power can be detected to 
find the .beginning of the frame. The carrier frequency offset 
is found in the frequency domain after applying a Hanning 
window and taking the FFT. The null symbol is also used in 
[M]. This extra overhead of using a null symbol is avoided by 
using the technique described in this paper. If instead of a con- 
tinuous transmission mode, a burst mode is used, it would be 
difficuh to use a null symbol since there would be no difference 
between the null symbol and the idle period between bursts. 

Van de Beek [3] describes a method of using a coirelaiion 
with the cyclic prefix to find the symbol timing. If this method 
were used to find the symbol timing, while using one of the 
previous methods to find the carrier frequency offset, there 
would still be a problem of finding the start of the frame to 
know where the training symbols are located. 

Classen introduces a method which jointly finds both the 
symbol timing and carrier frequency offset (5]. However, it 
is very computationally complex because it uses a trial and 
error method where the carrier frequency is incremented in 
small steps over the entire acquisition range until the correct 
carrier frequency is found. It is impractical to do the exhaustive 
search and go through a large amount of computation at each 
possible carrier frequency offset. 

This paper introduces some modifications of Classen's 
method which t)oih greatly simplify the computation necessary 
for synchronization and extend the range for the acquisition 

U(»()-677!S/»i75IO.W) O 1997 IEEE 



Abstraci^A rapid synchronization method is presented for 
an orthogonal frequenc>'*di vision multiplexing (OFDM) system 
using either a continuous transmission or a burst operation over 
a frequency-selective channel. The pre.scnce of a signal can be 
detected upon the receipt of Just one training sequence of two 
symtMls. The start of the frame and the beginning of the s>'mbol 
can be found, and carrier frequency offsets of many subchannels 
spacings can be corrected The algorithms operate near the 
Cramer-Rao lower bound for Che variance of the frequency 
oftset estimate, and the Inherent averaging over many subcarriers 
allows acquisition at very low signal to-noisc ratios (SNR's). 

lnd€x Terms^ Carrier frequency, orthogonal frequency- 
division multiplexing, symbol timing csfliiiaClon. 



I. IfrrRODUCTION 

IN AN orthogonal frequency-division multiplexing (OFDM) 
system, synchronization at the receiver is one important 
step that must be performed. Tbis paper describes a method 
to acquire synchronization for either a continuous stream, of 
data as in a broadcast application or for bursty data as in 
a wireless local area tietwork (WLAN). In both cases the 
receiver must continuously scan for incoming data, and rapid 
acquisition is needed. The ratio of the number of overhead bits 
for synchronization to the number of message bits must be kept 
to a minimum, and low-complexity algonihms arc needed. 

Synchronization of an OFDM signal requires finding the 
symbol timing and carrier frequency offset. Symbol timing 
for an OFDM signal is significantly different than for a single 
carrier signal since there is not an "eye opening'* where a 
best sampling time can be found. Rather there arc hundreds or 
thousands of samples per OFDM symbol since the number of 
samples necessary is proportional to the number of subcanriers. 
Finding the symbol timing for OFDM means finding an 
estimate of where the symbol starts. There is usually some 
tolerance for symbol timing errors when a cyclic prefix is 
used to extend the symbol. Synchronization of the carrier 
frequency at the receiver must be pcrfomied very accurately, 
or there will be loss of orthogonality betuecn the subsymbols. 
OFDM systems are very sensitive to carrier frequency offsets 
since they can only tolerate offsets which are a fraction of the 
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Fig. I. Block diagram of OFDM Itanimilter. 
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Fig. 2. Block diagram of OFDM receiver. 

of earner frequency offset. The method in this paper avoids 
the extra overhead of using a null symbol, while allowing 
a large acquisition range for the carrier frequency offset. By 
using one unique symbol which has a repetition within half a 
symbol period, this method can be used for bursts of data to 
find whether a burst is present and to find the start of the bursi. 

Acquisition is achieved in two separate steps through the 
use of a two-symbol training sequence, which will usually 
be placed at the start of the frame. First the symbol/frame 
timing is found by searching for a symbol in which the 
first half is identical to the second half in the time domain. 
Then the carrier frequency offset is partially corrected, and 
a correlation with a second symbol is performed to find 
the carrier frequency offset. 

II. OFDM Principles 
The OFDM signal is generated at baseband by taking the 
inverse fast Fourier transfonn (IFFT) of quadrature amplitude 
modulated (QAM) or phase-shift keyed (PSK) subsymbols 
Ck = flfc (Fig. I). In the figure^ the block P/S represents 
a pa^llcl-to^scrial converter. An OFDM syrnbol has a useful 
period T and preceding each symbol is a cyclic prefix of length 
Tg, which is longer than the channel impulse response so 
that there will be no intersymbol interference (ISl) [61. The 
frequencies of the complex exponentials are fk = k/T, and 
the useful part for 2N + 1 subcarriers is given by 

Ck exp(j27rA0. 0 < t < T. (I) 

Jk=-.V • 

The baseband signal is quadrature modulated, up-converted to 
the radio frequency (RF) and transmitted through the channel. 

At the receiver (Fig. 2), the signal is down-converted to an 
intennediate frequency (IF), and quadrature demodulated. The 
block S/P represents a serial-to-parallel converter. A carrier 
frequency offset of causes a phase rotation of 27r^ A /. 
If uncorrected this causes both a rotation of the constellation 
and a spread of the constellation points similar to additive 
white Gaussian noise (AWGN). A symbol-timing error will 
have little effect as long as all the samples taken are within 
the length of the cyclically-extended OFDM symbol. 

III. Estimation of Symbol Timing 

A. Symbol Timing Estimation Aigorithpn 

The symbol timing recovery relies on searching for a 
training symbol with two identical halves in the time do- 



main, which will remain identical after passing through die 
channel, except that there will be a phase difference be- 
tween them caused by the carrier frequency offset. The two 
halves of the training symbol are made identical (in time 
order) by transmitting a pscudonoise (PN) sequence on the 
even frequencies, while zeros are used on the odd frequen- 
cies. This means that at each even frequency one of the 
points of a QPSK constellation is transmitted. In order to 
maintain an approximately constant signal energy for each 
symbol the firequency components of this training symbol are 
multiplied by at the transminer. or the four points of 
the QPSK constellatioii are selected from a larger constel- 
lation, such as 64-QAM, so Chat points with higher energy 
can be used. Transmitted data will not be mistaken as the 
start of the frame since any actual data must contain odd 
frequencies. Note that an equivalent method of generating 
this training symbol is to use an IFFF of half the normal 
size to generate the time domain samples. The repetition is 
not generated using the IFFT, so instead of just using the 
even frequencies^ a PN sequence would be transmitted on 
all of the subcanriers to generate the time domain samples 
which are half a symbol in duration. These time-domain 
samples are repeated (and property scaled) to form the first 
training symbol. 

The second training symbol contains a PN sequence on the 
odd frequencies to measure these subchannels, and another 
PN sequence on the even frequencies to help detennine 
frequency offset. Table I illustrates the use of PN seqtimces 
in the training sequence for an OFDM signal with nine 
subcarriers with the points chosen from a subset of a 64- 
QAM constellation. The selection of a particular PN sequence 
should not have much effect on the performance of the 
synchronization algorithms. Instead the PN sequence can be 
chosen on the basis of being easy to implement or having a low 
peak-io-average power ratio so that there is little distortion in 
the transmitter amplifier. 

Complex samples are taken by mixing the received 
signal down to the IF, splitting the signal into two branches, 
multiplying by both the in-phase and quadrature local oscil- 
lators, and low-pass filtering and sampling to get baseband 
in-phase and quadrature components (Fig. 2). This can be 
e.\pres.sed mathematically by writing the IF local oscillator 
for the in-phase branch as 



(2) 
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TABLE I 

lU-eSTRATlON OF USE OF PN SEQL ENCES FOR TraIMNC SvmBOLS 



freq. num. k 




C2,k 




-4 








-3 


0 


-5-5j 




-2 


-7+7j 


-5-5j 


i 


-1 


0 




0 


7+7j 


-5-5j 


-I 


1 


0 


5+5j 




2 


7-7j 


--5+5j 




3 


0 


5-5j 




4 


7+7j 


5+5j 


1 



and the IF local oscillator for the quadrature branch at the 
receiver as 



- 2sin(23r/iF0- 



(3) 



Let the output of the mixer after down-conversion be s(t). The 
demodulated signal before the sampler can be expressed as 

r(0 = LPF{5(0 Vi{i)) + jLFF{s(t) v,(i)} (4) 

where LPF{-} means to low-pass filler the terms in the 
argument. The output of the in-phase branch is considered to 
be real and the output of the quadrature branch is considered to 
be imaginary. This is a mathematical convention to represent 
the in 'phase and quadrature components as a complex number. 
After sampling, the complex samples are denoted as rm- 

Consider the first training symbol where the first half is 
identical to the second half (in tinve order)« except for a phase 
shift caused by the carrier frequency offset. If the conjugate of 
a sample from the first half is multiplied by the corresponding 
sample from the second half (T/2 seconds later), the effect of 
the channel should cancel, and the result will have a phase 
of approximately ^ = ttTA/. At the start of the frame, 
the products of each of these pairs of samples will have 
approximately the same phase* so the magnitude of the sum 
will be a large value. 

Let there be L complex samples in one-half of the first 
training symbol (excluding the cyclic prefix), and let the sum 
of the pairs of products be 



which can be implemented with the iterative formula 
Pld^ 1) = Fid) 4. (rJ+tr^^2L) - (r>rf+t)- 



(5) 



(6) 



Note that d isn time index corresponding to the first sample in 
a window of 2L samples. This window slides along in time as 
the receiver searches for the first training symbol. The received 
energy for the second half-symbol is defined by 



<7) 




-0.5 0 0.5 1 
Timing offset (symbols) 

Fig. 3. Example of ihe liming metric for the AWQN channel (SNR = 10 
dB). . 

which can also be calculated iteratively. R(d) may be used 
as pan of an automatic gain control (AGC) loop. A timing 
metric can be defined as 



M{d) = 



\Pid)\\ 



(8) 



Fig, 3 shows an example of the timing metric as a window 
slides past coincidence for the AWGN channel for an OFDM 
signal with 1000 subcarriers, a carrier frequency offset of 12.4 
subcarrier spacings, and an signal-to-noise ratio (SNR) of 10 
dB, where the SNR is the total signal (all the subcarriers) to 
noise power ratio. The timing metric reaches a plateau which 
has a length equal to the length of the guard inter>-al minus 
the length of the channel impulse response since there is no 
ISI within this plateau to distort the signal. For the AWGN 
channel, there is a window with a length of the guard interval 
where the metric reaches a maximum^ and the start of the frame 
can be taken to be anywhere within this window without a loss 
in the received SNR. For the frequency selective channels, the 
length of the impulse response of the channel is shorter than 
the guard interval by design choice of the guard intervaU so 
the plateau in the maximum of the timing metric is shorter 
than for the AWGN channel. 

This plateau leads to some uncertainty as to the start of 
the frame. For the simulations in this paper, OFDM symbols 
are generated with 1000 frequencies, -500 to 499. They are 
slightly oversampled at a rate of 1024 samples for the useful 
part of each symbol. In an actual hardware implementation, the 
ratio of the sampling rate to the number of frequencies would 
be higher to ease filtering requirements. The guard interval is 
set to about 10% of the useful part, which is 102 samples. 

B. Performance of Symbol Timing Estimator 

There arc two issues to consider when evaluating the 
performance of the symbol timing estimator. First, since the 
timing metric is also used to determine whether the training 
sequence has been received, there is a probability of either 
missing a training sequence and not detecting the signal or 
falsely detecting a training sequence when none is. there. In 
this paper the distribuiion of the timing metric at the correct 
Stan of the frume is calculated. Using this distribution, the 
number of samples that need to be proces.sed during the 
detection phase can be determined and a threshold can be 
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set for this detection. Second, there is some degradation in 
perfonnance if the symbol timing estimate deviates from the 
correct region. Simulations are performed lo find the effect of 
extra interference that is introduced by poor symbol taming 
estimates for two types of channels. 

I) Distribution of Tuning Metric: Let each complex sam- 
ple = s,n + n^M be made up of a signal and a noise 
component. Let the variance of the real and imaginary com- 
ponents be: 



(9) 
(10) 

so that the SNR is al/al. To find the mean and variance of the 
estimator at the best symbol liming, first look at P(d) which 
can be written as 

. At Ihc coirect symbol timing, this can be broken into pans that 
are in-phase and quadrature to the s'^^^^^^sa^^+m^i^ product 
which has phase 0. This is just another way of looking at 
the problem with a new set of axes with one axis in the ^ 
direction and the other axis perpendicular to it. For usable 
values of SNR, when the magnitude is taken the quadrature 
part will be small compared to the in-phase part and can be 
neglected, so 

£-1 



(12) 

where inPhase^f } means the conij>oncni in the 0 direction. 
The quadrature part is neglected because the Rician distri- 
bution can be approximated by a Gau.<csian when taking the 
envelope of a dominant signal with Gaussian noise [8]. Tlie 
dominant in-phase part is so much larger than the quadrature 
part thai the envelope can be approximated by the in-phase 
part. Note that the first temi of \P(d^^^)\ is much larger 
than the second tenm since the ^:.^..^5u,,,+m+L produces 
all have phase 4> and add in-phase. v^hile all the other terms 
add with random phases. By the cenirjl limit theorem (CLT), 
|F(<fopt)| is Gau.vsian with mean 2Lff; since each of the L 



temis </,p.+,„^Jop»+m+£ has an expected value of 2a^ and 
all the other terms have an expected value of zero. 

For /?(dopt). the magnitude of each term could be taken by 
adding the squares of the real and imaginary parts. Instead, 
we can define a new set of orthogonal axes in which one axis 
is in the direction of the term 5 j,^, and the other is 
perpendicular to it. 

- H (l*'^p'+m+L|+inPhasc,,^^^^^^,{n,,.^, + ^4t})^ 

m:=0 



+ (quadratures^ Jn^^,,^^^^}) = 



(13) 

Here. i'^Phase,^^^^^^^^ {•} means lo take the component in 
the direction of . By the CLT. R{d) is also Gaussian 

with mean 2i(af -|- al). Note that for usable values of 
SNR, the mean is much greater than the standard deviation. 
Thus, for the Gaussian approximation, the probability of 
^(^opi) < 0 is insignificantly small making the approximation 
of the nonnegalive i2(^fopi) by a Gaussian reasonable. Another 
equivalent way of thinking about the distribution is that 
•'^(^^opt) is Rician with the mean much larger than the standard 
deviatiorl. In this case the Gaussian approximation may be 
used [8J. 

Define the square root of M{d) to be q{d) = \P{d)\/R(d) 
(see (14) at the bottom of the page). 

Since the standard deviations of both the numerator and 
denominator are much smaller fhan the means, the approx- 
imation aCl + b)/{l + c) % a(l + 6 - c) can be used. 
Another way to explain this Is that there is a Gaussian random 
variable in the numerator and a Gaussian random variable in 
the denominator, and since the standard deviation of both of 
these Gaussian random variables is much smaller than their 
miean values, tlie above appro.ximation can be used to write 
their ratio as a single Gaussian random variable. As long as 
this approximation holds, 9(c/»pi) is Gaussian with 

This can be justified because linear operations on a Gauss- 
ian random variable will result in another Gaussian random 
variable (9). 

When calculating the variance, note that 



and 



inPhasc^ ( s j^^^ ^ ^ n,/,^, ^^^l] 



*?(rfoni) = 



(14) 
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Since these tenms are the same in both the numerator and 
denominator, they do not contribute to the overall variance. 
Then 



(17) 



M[dotn) is {qid^pt))*, which is the square of a Gaussian 
random variable. Since the variance is much smaller than the 
mean, this also can be approximated by a Gaussian random 
variable since linear operations performed on a Gaussian ran* 
dom variable will produce another Gaussian random variable: 

M(d,pt)=(M,-t-n(0,a|))' 

= /i| + 2M,n(0,<T2)+ (r»{0,«T|))' 

«p| + 2p,n(0,<r|) (18) 

where n(^, cr^}. is used to denote a Gaussian random variable 
with a mean of p and a variance of o^. 
The expected value is 



IXM = E(A/((t,p.)l = 



(19) 



(20) 



and the variance is 

At high SNR the mean is approximately 1 and the variance is 
approximately 2aJ [(1 + l)a?a2]/[L(irJ)^] « 4/(Z • SNR). 

The value of A/((fopt) also can give an estimate of the SNR. 
which is 



1- VA'^(dopt)' 



(21) 



This equation is derived from (19). This can be written as 
SNR = -1 4- 1/(1 • 9(^f>i))- The denoniinaior (1 - v(^pt)) 
is Gaussian, and its reciprocal will also be Gaussian since 
the standard deviation is much smaller than the mean for low 
values of SNR. Again, this is a linear operation on a Gaussian 
random variable (after the approximation). Then 

^2 



£(SNR] = ^ = SNR 



var[SNRJ ^ 



(22) 



This estimator works well for the SNR below 20 dB. Above 
this level, Af (Jopt) is so close to 1 that an accurate estimate of 
the SNR can not be determined* but onjjMhal the SNR is high. 
For example, if M(d^pt) = 0.9, then SNR = 18.5 = 12.7 dB. 
This can be used to set a threshold so that very weak signals 
will not be decoded, or it can be used in a WLAN to feed back 
to the transmitter to indicate what data rate will be supported 
so that an appropriate constellation and code can be chosen. A 



lookup table can be implemented based on A/(r<„pt). so thai 
no square roots or divisions need to be performed. 

Even if there is a frequency selective channel, all the signal 
energy will go into the signal component temi except when 
the length of the channel impulse re.«;ponse becomes so large 
that it is longer than the cyclic prefix. At* this point, the energy 
located at longer delays becomes interference and would be 
added to the noise terms. 

At a position outside the first training symbol, the terms in 
the sum PidouiMe) add with random phases since there is not 
a periodicity for samples spaced by L samples. For the purpose 
of computing the mean and variance of the timing metric, the 
samples can be considered to be composed of just noi.se terms 
since noise terms will also add with random phases, so the 
statistics will be independent of the SNR. This assumption 
is verified by the simulation results shown in Fig. 6 which 
will be explained later in this section. The real and imaginary 
parts of P(douuidc) are Gaussian by the CLT. The sum of 
the square of two zero-mean Gaussian random variables, each 
with a variance of 1 is a chi-square random variable with 
two degrees of freedom and is represented by the symbol xl- 
The mean of xl ^ sui^ variance is 4 (9J. To simplify 
the computations, let the variance of the real and imaginary 
components of be: .. . 



(23) 



Both Re{P(dout»ide)} and Im{P(dout»w«)) arc Gaussian- 
distributed with zero mean and a variance of 2L<t^- Incor* 
porating this scaling factor, we have 



|P(rfout,ide)|^ = 2Latxl 

The mean and variance of |/*(dout«id«)P are 

£:[|P(<iout.idc)fl=4La; 
var(|P(£fo««idOli ^^L^^l 



(24) 



(25) 
(26) 



The denominator R{doutn\de) has a Gaussian distribution by 
the CLT. and its square is also Gaussian because the standard 
deviation is much smaller than the mean. Again, this is a 
linear operation on a Gaussian random variable (using the 
approximation), so the result is also Gaussian. Thus, 

/?(rfoatMcie)' - n{L^{2a^,)\L^2al)') (27) 

where the ^ operator means **is distributed.*' The ratio of these 
two random variables (after dividing both the numerator and 
denominator by a constant) is 



(28) 



Here, the variance of the Gaussian random variable is propor- 
tional to and can be neglected. The mean and variance 
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Fig. 4. £xpecied value of uming metric with I, = 3 12. Dashed lines indicate 
three standaid deviations. 
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Fig. 5. Mean and variance at the conrccc timing point with £ s 512 for Che 
exponential channel. 



of the timing metric are 



£?|Af (rfoufide)] = J 



(29) 
(30) 



var(A/(douu«ie)l = 

Fig. 4 shows a plot of the expected value of timing metric 
M(d) versus SNR at both the best liming instant and a point 
outside the training symbol. The dashed lines indicate three 
standard deviations from each curve. 

Figs. 5 and 6 show the results of simulations performed 
on the exponential frequency selective channels with 1(X)0 
subcarricrs for the correct timing and an incorrect timing, 
respectively. The conect timing point was chosen at the 
start of the useful part of the first training symbol, and the 
incorrect timing point was chosen one symbol after that. For 
each SNR value, the simulation was run 10000 times, each 
time generating different PN sequences, channels, and noise. 
Note that the simulated variance is slightly higher than that 
calculated by theory for the incorrect timing. This is caused 
because there Is a correlation between the signal terms r;^^,,, 
and while the calculation assumes that they are 

independent. 

2) Probability of Missing the Training Sequence or of False 
Detection: The probability distnbuiions calculated in Section 
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— ffiaery 







SNR(dB| 
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Fig. 6. Mean and variance ai an incorrect timing point with L =: 512 for 
the exponential channel. 



in-B-1 can be used to determine both the probability of not 
detecting a training sequence when one is present and of 
falsely detecting a training sequence when one is not present. 
As an example of how this can be done, let there be 1000 
subchannels and let L = 512. If the system is designed to 
detect a signal if the SNR is at least 10 dB, then from Eqs. 
(19) and (20). when the signal is present the mean is 0.827 
and the variance is 5,58 x 10'*, so the standard deviation is 
2.36 X 10"^ at the correct timing. If no signal is present, then 
from (28) the mean is 1.95 x 10*^ and the standard deviation 
is also 1.95 x 10"^. If the threshold is set at 0.1, then the 
margin for error when the training signal is present is 
0.827 - 0,1 ^ ■ 

2 36 x 10-^ standard deviations. (31) 

Similarly, when the signal is not present, the margin for error is 
0.1 - 1.95 X 10-"^ 



1.95 X 10-3 



= 50 standard deviations. 



(32) 



If the desired probability of error of missing a training se- 
quence is 1 X lO'^, this requires that the threshold be 3.1 
standard deviations below the mean found with (19). When 
computing the probability of false detection, note that the 
training signal is not present for most of each frame. If there 
are 100 OFDM symbols within one frame, then most of the 
time die training symbol is not present. Since the sliding 
windows for the symbol timing estimator are half a symbol 
long, there can be about 200 uncorrelated values of the symbol 
timing estimator within one frame. If the probability of false 
detection within one frame is desired to be I x 10'^, then the 
probability of false detection at any point in lime should be 
about 5 X 10"^, From [ 10] the cumulative distribution function 
for xi is 

^U5) = l-oxp[-.x5/>]. (33) 

To get a probability of error of 5 x lO"** requires ^ 
less than 24.4, which corresponds to 

(24.4 - 2)/2 = ! 1.2 standard deviaiion.s. (34) 
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Fig. 7. Relauofiship of signal, noise, and intetference power to the symbol 
liming position for the AWCN channel. The shaded portion in Che fiat plot 
indicates the range where the synchronizer can estimate the start of the symbol 
to maximize SNIR. 



In this example, both the probabilities of missing a training 
sequence or of false detection within one frame are less than 
1 X 10"^. This allows the threshold to be adjusted to operate 
at a lower SNR if that is desired. Another option is to reduce 
the amount of computation performed while searching for the 
training sequence by not processing every sample. This could 
be useful in a burst mode when data may not arrive vciy often. 

3) Reduction in SNIR: TWo channels are used in Simula- 
tions to measure the performance of the estimation algorithms. 
First the AAVGN channel is used to show how the algorithms 
perform on a benign channel, and it is also used as a basis 
for comparison. Then a hequency selective channel with an 
exponential power delay profile is used to show that the 
algorithms perform well in a more realistic environment. 
Sixteen paths are chosen with path delays of 0» 4» 8p* • 60 
samples, so the impulse response of the channel is shorter than 
the guard interval. The amplitude of each path is calculated 
from the exponential distribution, so 

i4i = exp{-r,/60) (35) 

where A; is the amplitude of the ith path and r,- is the delay 
of the ith path is samples. The phase of each path is chosen 
from a unifonn distribution from 0 to 2ir. 

The reduction in the signal-to-noise-plus-interference ratio 
(SNIR) due to a symbol timing error can be calculated for the 
simulated channels. SNIR can be defined as 

_2 



SNIR = 



(36) 



where a; is the variance of the intersymbol interference (ISI) 
added by incorrect symbol timing. Fig. 7 plots the average 
signal, noise, and interference power levels versus time for an 
AWGN channel. Note that interference here refers to the ISI 
from symbols before and afier the training symbol in the time 
domain. The time from -Tg to 0 is the length of the guard 
interval, and the time from 0 to T is the length of the useful 
part of the OFDM symbol. If the synchronizer estimates that 
the useful part of the symbol sturts at any time within the guard 
interval, which is the shaded area in the first plot of Fig. 7, 
then there is no reduction in SNIR due to incorrect symbol 
timing. However, if the synchronizer estimates that the start 
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of the symbol is outside the guard interval, there will be both 
a decrease in signal energy and an increase in interference, 
resulting in a lower SNIR. This occurs because samples from 
the previous or next symbol are input into the FFT along with 
samples from the current symbol. For example, if the start of 
the symbol is estimated to be at cither time T/11 or at time 
-T, r/11. then die signal-to-interfcrence ratio (SIR) is 10, 
assuming that the interference level is the same as the signal 
level, as would t>e the case with a continuous transmission 
but not with bursts of data. The SNIR will be lower than 
10 when the noise is included. For the simulations in this 
paper, it is assumed that a continuous transmission format is 
used because this represents the worst case in terms of self- 
interference. With multipath delay spread, the symbol timing 
tolerance window from -Tp to 0 is reduced by the length of 
the channel impulse response. 

Two methods to detennine the symbol linmg are coinpared 
on the basis of reduction in SNIR. The first method is to 
simply find the maximum of the timing metric. The second 
method is to find the maximum, find the points to the left and 
right in the tirne domain, which are 90% of the maximum, 
and average these two 90% times to find the symbol timing 
estimate. The rationale behind this method is that the best 
timing points typically He tn a plateau. By trying to determine 
the center of this plateau, it is more likely that the estimate 
will not fall slightly off the plateau. Table U compares the two 
methods on the basis of reduction in SNIR for the AWGN and 
frequency selective channel. For each type of channel. 10 OCX) 
simulations are run at each SNR, and in each run. different PN 
sequences, channels, and noise are generated. The averaging 
method performs significantly better than simply finding the 
peak of the timing metric, and it involves only slightly more 
computation. With a 40-dB SNR, a symbol timing off set of one 
sample away from the plateau would result in about 10 dB in 
degradation in SNIR for the AWGN channel. The averaging 
method of finding the symbol time resuhed in no degradation 
in 10000 runs for the AWGN channel and a degradation of just 
under 0.06 dB for the exponential channel at an SNR of 40 dB. 

IV. Estimation of Carrier Frequency Offset 

A. Carrier Frequency Offset Estimation Algorithm 

The main difference between the two halves of the first 
training symbol will be a phase difTerence of 

0 = ^T£iJ (37) 

which can be estimated by 

*> = auglo(P(d)) (38) 
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near the best timing poini. If |^| can be guaraniced to be less 
than jr, then the frequency offset estimate is 

A/ « 4,/{nT) (39) 
and the even PN frequencies on the second training symbol 
would not be needed. Otherwise, the actual frequency offset 



would be 



(40) 



where z is an integer. By panially correcting the frequency 
ofTsei^ adjacent carrier interference (ACI) can be avoided^ 
and then the remaining offset of 2z/T can be found After 
the two training symbols are frequency corrected by ^/(nT) 
(by multiplying the samples by exp(-j2i<^/T)), let their 
FFT s be xi^k and X2,k» and let the differentially-modulated 
PN sequence on the even frequencies of the second training 
symbol be vjt (a^ illustrated in Table I). The PN sequence 
Vk will appear at the ou^jut except it will be shifted by 2z 
positions because of the uncompeiisated frequency shift of 
2z/T, Note that because there is a guard interval and there 
is still a frequency offset, eveii if there were no differential 
modulation between training symbols 1 and 2 (i.e., PN se- 
quence Vk = 1), there would still be a phase shift between 
zi^k and X2,* of 27r(r + Tg)2z/T. Since at this point the 
integer z Is unknown, this additional phase shift is unknown. 
However, since the phase shift is the same for each pair 
of frequencies^ a metric similar to (8) can be used. Let X 
be the set of indices for the even frequency components 
X = {^W,^W + 2,....~4,^2,2,4,..«,lV--2,W^}.Thc 
number of even positions shifted can be calculated by finding 
g to maximize 

with integer y spanning the range of possible frequency offsets 
and W being the number of even frequencies with the PN 
sequence. Then the frequency offset estimate would be 



(42) 



B. Performance of Carrier Frequency Offset Estimaior 

Since the carrier frequency offset is found in two steps, by 
finding the fractional and then integer pan of the offset, there 
are two issues to consider in evaluating the carrier frequency 
offset estimator. First the variance of the estimator of the 
fractional part of the frequency offset can be found. Second, 
the distribution of the carrier offset metric can be found for 
both the correct and incorrect integer offsets to determine the 
probability of making an error in finding the integer part of 
the carrier frequency offset. 

/; Variance of Fractional Part of Carrier Frequency Offset 
Estimator: Using the method in (2] 

Since the carrier frequency offset estimate (in subcarrier 
spacings) is made up of the sum of the initial e.stimate and 
an even integer, the variance of the initial estimate. />/ir. will 
be equal to the variance of the final e.<timate if the integer g is 
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Rg. 8. Comparison of Moose's method and neu- method of esiimaiine 
frequency offset. 

equal to ^cprrect- The effect of the remaining frequency offset 
on the system performance is discussed in (2] and will not be 
discussed here. By using the method in (7J, the Cramer-Rao 
bound is found to be 

This is not surprising because Moose shows that his estimator 
is the maximum-likelihood estimator (MLE) of differential 
phase, and Rife states that the Cramer-Rao bounds arc almost 
met by the MLE with high SNR. Since the frequency offset 
estimate is made by averaging over hundreds or thousands of 
subcarriers, the effective SNR is usually very high- To illustrate 
that the frequency acquisition range is greatly widened with 
thU new ihethod, Fig. 8 compares the error variance for 
Moose's frequency estimation methods (with two repeated 
half-symbols) with the new method for 1000 earners with an 
SNR of 10 dB. The simulations were performed with lOOOO. 
runs per firequency offset value. Since Moose^s method is 
designed to work only with a very small frequency offset, 
it fails for larger frequency offsets. In Fig. 8, the Cramer-Rao 
bound is 2.03 x IQ-^. Note that the simulation here for 
Moose*s algorithm uses the shortened symbols as described 
in f2J, which results is a frequency acquisition range of il 
subchannel spacing. 

Fig. 9 compares the error variance for the carrier frequency 
offset estimate to the Cramer-Rao bound for different values 
of SNR. In the simulation the carrier frequency offset was set 
to 2.4 subcarrier spacings so the estimator has to get the integer 
part correct and gel a good estimate of the fractional part of 
the carrier frequency offset to make a good overall estimate. 
Again, 1 0000 trials were run at each SNR value, 

2) Distribution of Frequency Offset Metric: The expected 
values and variances of the estimator B(g) can be calculated 
tn a similar method as in Section 111-8-2 lo determine if 
Orowract can be found reliably. At the correct frequency offset, 
all the signal products 5l.,.^.,_.u;s2.h>i,,.,_. have the 
same phase and 
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Fig. 9. Comparison of carrier frequency offset estimaie to Cramer-Rao 
bound. 
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Rg. 10. Expected value of carrier frequency ofTset metric with IP = 500 
Dashed lines indicaie three standard deviations. 

At high SNR the mean is approximately 1, and the variance 
is apprgximalely 

At an incon-ect frequency offset the signal products no 
longer add in phase, and B(^uicorrect) has a chi-square dis- 
tribution with two degrees of freedom with 

^(^(..n.o„«01=2^(l + ;j^)<i, (47) 

1 + ^fLtifld^ ^ _L 

(48) 

Fig. 10 shows a plot of the expececd value of frequency 
offset metric D(f)) versus SNR for both a correct and incorrect 
frequency offset for the 1000 subcarrier OFDM signal. The 
dashed lines indicate three standard deviations from each 
curve. 
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upon the receipt of just one training sequence. By avcraeine 
over al! the subchannels, it works wcU in frequency scleciive 
channels. This method also gives very accurate estimates of 
symbol timing and carrier frequency ofTset and provides a 
very wide acquisition range for the carrier frequency offset 
h also provides an estimate of tlie SNR, and the probability of 
false locks or missing the Uaining symbols is very low. For a 
wireless LAN, siich a fast and low-ovcrhcad synchronization 
process is necessary because there will be only one training 
sequence transmitted in each burst for synchronization. 
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V. Conclusion 



A method has been presented for the rapid and robust ia ^ ^ . . 

synchronization of OFDM signals, and acquisition is obtained ^'^^l^^^-^i,!:.:^:^ 
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